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Abstract 

Reheating after R 2 inflation proceeds through gravitational particle production of 
conformally noninvariant fields. We argue that the nonvanishing expectation value of 
flat directions generic in supersymmetric theories break conformal invariance of the 
fields coupled to them in a position-dependent manner due to quantum fluctuations. 
As a result modulated reheating can occur after the supergravity i? 2 -inflation. The 
resultant curvature fluctuation is a mixture of the one produced during inflation and 
that produced by modulated reheating. The spectral index takes a value between 
n s = 0.960 and 0.983, the nonlinearity parameter of the local-type non-Gaussianity 
can be f^L ~ ±10, and the tensor-to-scalar ratio is r < 4 x 10~ 3 . 

RESCEU-5/13 
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1 Introduction 



The R 2 inflation [lj is a unique inflation model in the sense that it is one of the 
oldest models of inflation [2] but still observationally viable. In fact, it is not only 
viable but preferred by the observational results of WMAP9 |3j with its predictions 
of scalar spectral indexjl n s = 0.964, and the tensor-to-scalar ratio, r = 3.9 x 1CP 3 , 
occupying the center of the likelihood contour. In order to further probe or falsify 
this model, it is desired to have more predictions that can be tested soon. Among 
such quantities is the nonlinearity parameter l of curvature fluctuations [5j which 
measures the deviation of their statistical distribution from Gaussian. 

Previously, R 2 inflation has been known to yield highly Gaussian curvature fluc- 
tuations, which can be easily understood regarding the extra scalar degree of freedom 
in /(i?)-type theory, dubbed scalaron <p, as the driving force of inflation [6]. Indeed 
(p has a very flat effective potential ideal for chaotic inflation [7] with a canonical 
kinetic term. In this paper we argue that nontrivial processes in the reheating phase 
after inflation may produce additional curvature perturbations, which can be appre- 
ciably non-Gaussian, through the so-called modulated reheating scenario [8j known 
to produce Jnl of local type [9], [10] . 

In this theory, inflaton is followed by damped oscillation of the scalar curvature 
R(t), which induces gravitational particle production of conformally non-invariant 
fields thereby reheating the Universe. This process is equivalently described by 
the decay of the scalaron-inflaton ip which oscillates around the potential minimum 
after inflation. In the scalaron picture, <p is coupled with conformally non-invariant 
fields only, so it is coupled with neither massless fermions nor gauge bosons (except 
gauge conformal anomaly [llj). Therefore, the previous study of reheating after R 2 
inflation [121 fl3"l [14"] including more general f(R) inflation [15] mostly considered 
creation of massless minimally coupled scalar bosons. 

However, since mass terms break the conformal invariance, scalarons can decay 
into massive particles through their mass terms |16] (see also \17\ [T8] for earlier 
works). These channels open up a new possibility to produce curvature perturbations 
because the mass terms, which are determined by the Higgs field, turn out to be 
spatially dependent reflecting quantum fluctuations of the Higgs condensation [19] 
generated during inflation. As a result the modulated reheating [8] may take place 
which contributes to the observed curvature perturbations. 

Unfortunately, however, such a mechanism does not work in the original R 2 
inflation model which is realized by adding a sing le term, i? 2 /(6M 2 ), to the grav- 
itational Lagrangian. In this simplest model the scalaron mass is determined (or 
bounded from above, in case there is other sources of curvature perturbation besides 
the scalaron's fluctuations) by the observed amplitude of curvature perturbations as 
M ~ 3 X 10 13 GeV. As a result, its decay rate per mode Ti = M 3 /(192vrA^ 1 ) ~ 8 
GeV is so small that the Universe would not be reheated until long after inflation. 

1 The latest observational results of Planck [32, 33 also support the R 2 inflation, with tighter constraints 
on n s , r and /at l ■ We discuss their results further in subsection 14.11 . 

2 The first prediction of the scalar spectral index for R 2 inflation was made correctly by Mukhanov and 
Chibisov [4] 
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During this long period of scalaron oscillation, the Higgs condensation also oscillates 
and decreases its amplitude. Therefore, by the time the reheating occurs, the stan- 
dard model (SM) Higgs field value would be too small to leave any observational 
trace in contrast to the case studied in [22]. So much is the story of R 2 inflation in 




On the other hand, it is barely likely that the SM is the ultimate theory with 
which we should describe the birth and evolution of the Universe. As the primary 
extension of the SM, there are a number of sound motivations to introduce the su- 
persymmetry (SUSY). Then the R 2 inflation - or its analogue - should be realized 
in the context of local supersymmetry or the supergravity (SUGRA). An attempt to 
realize R 2 inflation in this context was first made by Ketov [20J , and was improved by 
Ketov and Starobinsky [21]. Its most wonderful feature is that the scalaron mass in 
the reheating phase is not directly related to the amplitude of curvature fluctuations 
and can be much heavier than in the original model, so that the Universe can be 
reheated immediately after inflation by gravitational particle production. The other 
plausible feature is the existence of flat-direction scalar fields which have potentials 
much flatter than the SM Higgs field has, even after the effects of SUSY breaking 
and non-renormalizable superpotential are taken into account. Thus quantum fluc- 
tuations around a large expectation value of the flat direction can induce modulated 
reheating [23J, and generation of appreciable amount of local-type non-Gaussianity 
is possible in this model. 

The rest of the paper is organized as follows. In Sec.[2]we introduce R 2 inflation in 
SUGRA proposed in [21] and further analyzed in [23]. Then in Sec. [3] we incorporate 
a SUSY flat direction </> to analyze non-Gaussianity produced in this model. Sec. H]is 
devoted to conclusion. In the appendix we describe scalaron interactions with other 
fields and calculate its decay widths as an equivalent way to calculate gravitational 
particle production rate. 



Here first we briefly introduce higher curvature inflation model proposed by Ketov 
|20j and Ketov and Starobinsky [21] in SUGRA. The theory is defined by the action 




the SM. 



2 R 2 inflation in SUGRA 




(i) 



where F(JZ) is a function of the scalar curvature superfield 



Tl = ~B* - 6 (a a a b ^ ab - ia a ^B* + i^ a b a 



) 




(2) 



and £ is the chiral superspace density in a Wess-Zumino type gauge 



£ = e(x) [1 - 2i6a a i) a (x) + 96k 2 B(x)] . 



(3) 
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Here, e = y/^g, 4> a is the gravitino, B is an auxiliary scalar field, and k is the 
reciprocal of the reduced Planck scale Mpi = k~ 1 = mpi/\f&it = 2.4 x 10 18 GeV. 
Ignoring fermions, we find 



F(K) = F + F> (±R - ^BB* 



d 4 x e 



X S --B. 



-3XF {X*) + F' (X*) ( -R - 4XX* 



(4) 

+ h.c, (5) 
(6) 



Ignoring a pseudo-scalar partner (axion) of the scalaron, namely, taking X = X* 
and F = F*, the Lagrangian of the gravity sector reads, 



C = -6XF(X) + 2 Qi? - 4X 2 ^j F'(X), 



(7) 



together with the constraint which is nothing but the equation of motion (EOM) for 
the auxiliary field 



= 3F(X) + 11XF'(X) - [-R- AX 2 F"(X). 



1 



The choice of Ketov and Starobinsky is 



F(K) = ^f 1 K+^f 2 K 2 + ^f 3 K 3 , 



which leads the Lagrangian 



C = -hR + -f 2 RX + -f 3 R - 7h X 2 - 11/ 2 X 3 - 5/ 3 ^ 



and the constraint 



n - X 3 + ^X 2 -±(R E'u ) A' 



20/ 3 " 



30/ 3 



R, 



with 



i?n = 



21/i 

h 



(8) 



(9) 



(10) 



(11) 



(12) 



where f±, f 2 , and f 3 are positive constants. Note that we use the metric signature 
(— , +, +, +) so that our sign convention of 1Z and R = 12H 2 + 6H is different from 
that of Ketov and Starobinsky. 

Since the gravitational constant must be positive for stability of the system, 
F'(X) > and thus /1/3 > f 2 . The absence of ghost and tachyonic degrees of 
freedom requires f\ > and f 3 > 0. The sub-Planckian curvature during inflation 
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requires /3 > 1. The sub-Planckian scalaron mass after inflation requires /f S> fi 
in order to avoid large quantum gravity loop corrections. In summary, [21] 



h » l, hh > /I » h > 0. 



(13) 



In terms of the reduced Planck scale Mpi, the scalaron mass during inflation, 
M, and that in the reheating stage around the potential minimum, m, they can be 
expressed as 



fi = -M 2 Pl , h 



[m m pi 

2""- r " " vIV : 

respectively. Then the conditions (fT3|) read 



and h 



15M|; 

Hkf 2 "' 



m > 



1m, 

20 



M < M PU m <C Mp^ 



(14) 



(15) 



respectively 

In the high curvature regime R S> Rq 
reads 



2I/1//3 = 21M 2 /10, the Lagrangian 



+ J_ hR 2 + /2jR 3/2 
3 180 100 



6M 2 



+ 



3VW5 R 3 / 2 

Too 



■m 



(16) 



The dynamics of inflation realized in this Lagrangian has been studied by Ketov 
and Tsujikawa [23] in detail. They solved the EOM for the Hubble parameter of the 
Friedmann Universe, H(t), and expressed the number of e- folds, N, acquired during 
inflation after the epoch H(t) = H in terms of the dimensionless quantities 



a 



M 2 
mH 



and /3 



M 2 

w 



as 



N 



1 



126a 2 
8000/3 



3a (80^ - 21a - v / 7(63a 2 + 16000/3/3) 



4000/3 



In I ^ 63 " 2 + 800/3 ) + 21ay / 63a 2 + 1600/3^ 



+ 3 y 210 + 2^/35/3 

The square amplitude of curvature perturbation is given by 

100/3 



1250 ( MY 



3V35a + 



(17) 

1 In 2 + 3 In 5) 
(18) 

(19) 



We assume this curvature perturbation generated by inflation is responsible for the 
fraction A 2 of the total amplitude, V bs — 2.4 x 10 -9 , determined by COBE and 
WMAP [3] on the pivot scale K = 0.002MPC" 1 , namely, 



2.4 x 10~ 9 A 2 , 



(20) 
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where N* is the number of e-folds of inflation after the comoving pivot scale left the 
Hubble radius. 

The scalar spectral index n s and the tensor-to-scalar ratio r can also be expressed 
by a and j3 as 



n s = l — a-—, and r = (9V35a + 100/3V 4 , (21) 

100 3 2500 v H> v 1 



1 t^-oi and r = (9\/35i 

for the case A 2 = 1 [24"] . 

As is seen in (|16p . the original i? 2 -inflation is obtained in the small a limit, when 
N approaches 3//3 - 1/2 or /3 = 6/(2JV + l). Then the COBE-WMAP normalization 
yields 

= 7Mxl0 ~\ Mpi = 1.36 x 10- s Xv-lM Pl , v* = - + 1/2 (22) 
A* + 1/2 * 55.5 v ; 

which in turn determines the Hubble parameter as a function of N as 

ff= (* + l/2yV (23) 

Thus the Hubble parameter when the pivot scale left the Hubble radius is given by 

H* = 4.1 x 10~ 5 AiaT 1/2 Mpi. (24) 
On the other hand, the mass of the scalaron in the reheating phase is determined 

as 

M 2 a I/O i 

m = — — = 4.47 x 10" 6 A^ 6, a~ l M Ph (25) 

so that the scalaron mass in the reheating phase is unrelated with the CMB obser- 
vation, and can be chosen arbitrarily under the condition m <C Mp\. Here a* is the 
value of a at H = H*. Thus the decay rate of the scalaron (see Appendix A and set 
M s = m there), 

where A/5 stands for the number of scalar decay modes due to scalar kinetic inter- 
action terms, can be larger than the cosmic expansion rate at the end of inflation, 
H f ^ M/V12, if 



> 0.133A/" 5 1/3 vl 1/3 \ 1/3 Mpi = m rh . (27) 



In such cases, the Universe is reheated efficiently soon after inflation even if the 
inflaton-scalaron decays only through gravitational interaction. In what follows we 
fix vi = m r /, so that the Universe is reheated rapidly after inflation. 
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The above inequality ([27|) imposes an upper bound on a as 



a, < 1.6 x 10- 4 A 2 /V 7/6 (A/>/10 2 )V3. (28) 

Thus the rapid reheating after inflation automatically guarantees the original R 2 - 
inflation in the supergravity context. The scalar spectral index n s R2 and the tensor- 
to-scalar ratio rji2 from inflation are given by 

n sR2 = 1 - = 0.964, r m = = 3-9 X KT 3 , (29) 

respectively, where we have quoted the numerical values at N = 55. 

Finally we comment on the thermal history after inflation in this model. It 
has been known that finite amount of radiation is created after the onset of field 
oscillation even if only perturbative decay with the rate (|26p operates, and its energy 
density is given by p r (t) = |riI(i)Mp Z (see e.g. [27]). Since F is as high as Hj here 
by assumption, the radiation temperature just after inflation is comparable — in fact 
only 20% lower — to the case all the energy density is immediately converted to 
radiation at the end of inflation. Thus even if the preheating may occur just after 
inflation as analyzed by Ketov and Tsujikawa [24], the thermal history does not 
change significantly. Whether preheating occurs or not in the presence of large 
r deserves further study, though. The high temperature effect may modify the 
potential of the flat direction to initiate field oscillation due to thermal effects, which 
does not change 8(f) /(j) because both the homogeneous mode (f> and superhorizon 
long-wave fluctuations 5(f) would oscillate in the same manner. Furthermore, in our 
model, the scalaron decays rapidly before these oscillations dissipate field amplitude 
appreciably. So we may use the field values at the end of inflation below. 

Since the reheat temperature is quite high, the gravitino problem is a problem 
[25] , We should either adopt a model with large enough gravitino mass without 
changing cosmic expansion history, or invoke non-standard cosmology like thermal 
inflation [26]. In the former case, the pivot scale is pushed up to iV* = 60 due to the 
earlier radiation domination, while in the latter case it may be shifted to a value as 
small as iV* ~ 50. 



3 Modulated reheating through SUSY flat di- 
rections 

We now consider behaviors of a SUSY flat direction field in this inflation model. 
Although our favorite one is H u Hd flat direction in close analogue to the Higgs 
condensation in the SM, other combinations such as LH U or those involving three 
fields would also work with their appropriate choice of coupling constants, provided 
they do not suffer from Q-ball instability before the scalaron decay. 

Here we denote a generic flat direction by <f> and study its cosmological effectsU 

3 Note that a flat direction <fi is not the inflaton, but causes curvature perturbations by modulated 
reheating. The H u Hd flat direction can be associated with it by t H u oc (0, (f>) and t Hd oc (<fr, 0). 
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By definition (ft acquires a potential only through SUSY breaking and possible non- 
renormalizable terms in the superpotential, W D ^Jn-a with <]? = -^e lS ■ Since we 

are interested only in the (position-dependent) amplitude of the flat direction, we 
suppress possible CP violating A-terms to write the potential as 

1 <0n-l) 

^ 1 = 2 m ^ + ^M' n " 4 ' (30) 
M x 

for (ft <C Mpi. Here mo is the SUSY breaking mass assumed to be O(lTeV) or so, 
and Mx is a cutoff scale to characterize nonrenormalizable interactions. 

1 n-3 

The second term is dominant for (ft > m^~ 2 M^~ 2 = (j)NR which we assume is 
satisfied during inflation in the subsequent analysis. For n = 4 and 6 we find 4>nr = 
2.0 x 10- 8 Mpi and 1.4 x lCT 4 M Pi , respectively, for m = ITeV and M x = M Pi . 
With no Hubble induced mass terms here during inflation, the initial value of (ft 
is randomly distributed, and it does not evolve significantly during inflation if its 
effective mass is smaller than the Hubble parameter, V"[(/>] < Hj. This yields an 
upper bound on (ft as 



< 



(2n-2)(2n-3) 



(31) 



If it was not satisfied in the beginning, the field would roll down the potential until 
it is satisfied. Using (f2"2"j) and (f2"3"j) . we find 

1/2 



J.5 x 10- 4 AVV 1/2 Qg) M Ph (n = 4) 



(ft max = 2.5 x lO-^VVr 174 M PU (n = 6), (32) 

respectively. 

It is convenient to express solutions of the homogeneous mode, (ft, and fluctuation 
around it, Scft, as a function of the number of e-folds until the end of inflation, N, 



which is a decreasing function of cosmic time, using the Hubble parameter (|23|) [10 . 
The slow-roll EOM of (ft reads 

d$_ = 1 U = 6^ H ' = ^L (w) 

dN 3-H'/HH 2 (3iV + l)M 2 ' dN f 1 ' 

which is solved as 



cft{N) = eft, 



8(n-l)(n-2)^"- 4 ( N* + 1/3 

X 



M 2 M 2n-6 '"\N+ 1/3 



(34) 



Here (ft* is the initial value at N = N t ., which we take the epoch of the wave number 
of our interest, namely, the pivot scale A;* left the Hubble radius during inflation, 
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and a prime denotes differentiation with respect to N. In particular, at the end of 
inflation corresponding to N = 0, we find 



0(0) 
A 



'n- l)(n-2)^- 4 , , Ar 
1 + : i\ ^ « In (3iV* + 1) 



M 2 M] c n ~ 6 



(n - l){n - 2)4> 
M 2 M] c n ~ 6 



2n-4 



ln(3iV* + 1) < 



?? 



3(2n - 3) 



0/ =^(l + A)2n-4, 



In (3iV* + 1) = A max , (35) 



, max - 0.682 and A^ 



where we have used {3TJ and fl) ^ M/VU. We find A™=* 
0.757 for A 7 * = 55, respectively. 

Similarly, the super-horizon evolution of fluctuation can be obtained by pertur- 
batively solving the EOM 



1 



dN 



to yield 



6(fi(N) = 5cfi* exp 



x exp 



N 



H'/H 



1 



H 2 



+ 



Vn 



2H 2 



(36) 



, 3- H'/H H 2 

V dN' 
jv. 3 - H'/H H 2 



dN 



N 



1 



V. 



3 - H'/iT # 2 



exp 



■dJV 



N, 



3 - H'/i? i7 2 



v u dN n 



(37) 



Here the initial condition is given by 



- H Jf /{2 / K). At the end of inflation, we find 

2n-3 



where is given by 



e 



2n - 3 



,2n-4 



/2n-3 

(2n - 3)A 



,2n-4 
6 / 



2 



(38) 



2n-3 



(1 + A) 2+ ^^ 0/ 



(39) 



To calculate the power spectrum of curvature perturbation generated by modulated 
reheating, only the first term in (|38p is important, whereas the second term is also 
necessary in principle to calculate non-Gaussianity. 

Since the Universe is reheated rapidly after inflation, we can estimate the am- 
plitude of curvature perturbation produced by the modulated reheating, ( m , and 
its non-Gaussianity using 4>f and 5(f) t. For this purpose it is convenient to use the 
scalaron picture and use the formulas given in |10j . 
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First Cm is given by 



where Q(x) is denned by 



xQ'(x) 6 -^, (40) 



Q(x) = Q(T/H) = - A \n 



dN' 



-e 



47V' P<p( N 



H{N>) p vf 



(41) 



Here p^ is the effective energy density of the scalaron-inflaton and the subscript / 
denotes the value at the end of inflation [10J. Our choice of the mass parameter 
corresponds to x = 1, where we find Q'(l) = —0.0890 and Q"(l) = 0.126. Here the 
total decay rate T is given by T = Ms^sk (A/s = M a + A/#+ all other light scalars) 
and the space-dependent part of the rate is expressed as 



(Af a + Afip)y 2 m(j)5(f) _ N FL y 2 mc t 



2AttM 2 pi " 2AttM 2 pi 



(42) 



where we have put h = y and replaced M s by m in the formulas (|76p and (|77|) . Thus 
Cm is given by 

- Cm = 0.712 MFL f^ = 2.33 x lO^^X^^l, 
J\f s ' m 2 

-(ir(#)(i)" l/3 

Here p is a numerical factor normalized by typical values of the number of decay 
modes and Yukawa coupling. Note that the number of decay modes for a pair of a 
single species of quark and antiquark, say t and t, is = 12 and supersymmetry 
doubles it by the same contribution of Tskm- We also note that Yukawa coupling is 
typically large, which might even saturate perturbative bound y 2 /(4ir) ~ 1 at high 
energy scale of our concern [28]. Using ([35]) . (j38j) and (|24l) . we find 



- Cm = 1.53 x 10-V* 1/6 A 1/3 (1 + A)-^-^. (44) 

Mpi 

It should contribute to the fraction 1 — A 2 of the observed power spectrum by as- 
sumption, so V(- m = (1 — X 2 )V bs- Then 4>*/Mpi can be expressed as 

3.2 x lO-VV^A- 1 ^! _ A 2 ) 1 /2 (1 + A )*E§. (45) 



Mpi 



The spectral index, n sm , of perturbation produced by modulated reheating is given 
by 

nsm = 1 - = 0.982, (46) 
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where the last equality is for the case N = 55. 

Next we turn to the non-Gaussianity. It is known that modulated reheating 
produces an appreciable amount of non-Gaussianity of local type [91 [10], and the 
nonlinearity parameter f NL of Cm, fNLm, is given by 

fNLm ~ 6Q^) + + 3 " 7( " T) w^f (47) 

for x = 1. The full nonlinearity parameter is given by 

V L 2 2 

fNL = fNLmTffir = fNLm{^ ~ A ) 

ofes 

= 13.2(1 - A 2 ) 2 + 35.0(T - 1)(1 + A)~^=r^ 1/3 A 4 / 3 (1 - A 2 ), (48) 

where we have used ([2"7) . ([35]) and (gU). 
From dSSJ), (05} and (J22J) A satisfies 

A = 1.38 x 1Q 2 ln(37V +1) ^.4 + A) 6 A -io/3 (1 _ A 2 )2 ( MpiY (4g) 
In 166 \Mx / 

for n = 4, and 

A = 4.78 X 1Q -8 ^(3^ + 1) ^2/3 8 (1 + A) 10 A -14/3 (1 _ A 2 } 4 ( Mpi\ " (g(J) 

In 166 \ Mj / 

for n = 6. Here note that we are not allowed to take A = exactly, because we have 
first fixed M from (1221) and this procedure is invalid there. Since this problem arises 
if and only if we take A exactly equal to zero from the beginning, our formula for 
Jnl and n s below are applicable even in the limit A is small and approaches to zero. 

For definiteness, let us take z/* = 1 (i.e. N* = 55) and \i = 1 below. Then we can 
express /jvl as a function of A and A for each n as 

f NL = 13.2(1 - A 2 ) 2 + 35.0K(n, A)A 4 / 3 (1 - A 2 ), (51) 

where 

2n-3 4n-6 2n-3 

K(n, A) = (T — 1)(1 + A)~— = (2n - 3)A(1 + A)~^=* - (1 + A)~^~. (52) 

We can also determine the cut off scale from (j49j) and (jSTij) as 

M x = 12A" 1 / 2 (1 + A) 3 A~ 5/3 (1 - A 2 )M Pi , (n = 4), (53) 

M x = 6.0 x W^A-^il + A) 5 / 3 A- 7/9 (l - A 2 ) 2/3 M P ,, (n = 6), (54) 

respectively. 

Now let us evaluate K(n, A) for n = 4 and n = 6. We find 

- 1 < K(4, A) < -1.93 x 10~ 2 , (55) 
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which is monotonically increasing as A, and the inequalities are saturated at A = 
and A max = 0.682, respectively. For n = 6 we find 

- 1 < A'(6,A) < 0.327, (56) 

where the minimum is realized again at A = and the maximum at A = 0.451. 

The case A — > corresponds to the limit Mx is large and nonrenormalizable 
interaction is absent. Then, independent of n we find 

f NL = 13.2(1 - A 2 ) 2 - 35.0A 4/3 (1 - A 2 ), (57) 

which takes maximum /jvl = 13.2 in the limit A — > 0, and minimum Jnl = —7.9 at 
A = 0.749. The former corresponds to the case all the observed fluctuations are due 
to the modulated reheating with the spectral index n s = 0.982, while at A = 0.749 we 
find n s = 0.972 for = 1, because the spectral index of total curvature perturbation, 
n s is given by 

2A 2 + 2 

n s = X 2 n sR2 + (1 - A 2 )n sm = 1 — , (58) 

in general. 

On the other hand, with K(4, 0.682) = -1.93 x 10~ 2 , one finds 

f NL = 13.2(1 - A 2 ) 2 - 0.675A 4/3 (1 - A 2 ). (59) 

In this case /jvx decreases monotonically from 13.2 to —8.32 x 10 -3 as we increase 
A from to 0.988. Then f^L increases to as A — >• 1. 
For the case A(6, 0.451) = 0.327 we have 

f NL = 13.2(1 - A 2 ) 2 + 11.5A 4/3 (1 - A 2 ), (60) 

which takes maximum /tvl = 13.5 at A = 0.148 and then decreases to as A is 
increased. 

Our results suggest that as far as fj^L parameter is concerned, the nonrenor- 
malizable term in the potential of the flat direction does not have observationally 
significant effects, which means that our result is robust in this respect. 

If we take smaller values for fi, then the second term in (|48p contributes less 
and the overall /ni can become larger for n = 4 and smaller for n = 6 than the 
demonstrated case above. 

4 Conclusion 

In the present paper we have reconsidered cosmic history after i? 2 -inflation in SUGRA 
in which reheating proceeds through gravitational particle production of conformally 
noninvariant fields. We have argued that conformal invariance is broken through a 
nonvanishing expectation value of a SUSY flat direction which gives a position de- 
pendent mass for those fields coupled to it due to its long-wave quantum fluctuations 
acquired during inflation. As a result such fluctuations induce modulated reheating. 
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We have shown that with reasonable values of the model parameters, on condition 
that the universe is reheated rapidly, curvature perturbation produced by modulated 
reheating may contribute to the observed fluctuation significantly and yields a local 
non-Gaussianity with f^L ~ ±10. 

The intriguing feature of our model is that the observable parameters n s , r, /jvx 
are mutually correlated with each other through the ratio of the two contributions 
A 2 , although the exact one-to-one correspondence is difficult due to the uncertainties 
in the pivot scale and the particle-physics parameter ji. For example, from (|58p 
and 

2 48A 2 48A 2 

r = rR2X = (2Ar, + i )2 = mV' () 

one can obtain A and z/* (or N*), which can be compared with /jvi- In future, 
information on the thermal history carried by v„ can also be tested by space based 
laser interferometer such as DECIGO [29, 30J, and can be compared with our result. 



4.1 Note added after Planck 2013 results 

After we submitted the original version of this article to the arXiv, the Planck 
collaboration released their data. It gives [32j [33] 

n s = 0.9653 ± 0.0069 , r < 0.13 (95% CL), f NL = 2.7±5.8 (68% CL), (62) 

where we have quoted values from Planck combined with WMAP polarization and 
lensing data sets for n s and r. Their results are basically in good agreement with 
the predictions of our model. 

To be precise, the lower value of n s , which favors larger A, does not yield any 
sizable non-Gaussianity with our model. On the other hand, the higher allowed 
value like n s = 0.967 yields 

r = 3.2 x 10~ 3 , f NL = {-4.8, 0.28, 2.1} (63) 

for {K(n,Q + ),K(4,QM),K(6,0A5)}, A = 0.91, p = 1 and u m = 1. Thus our model 
covers lower allowed range of the observed /nl, namely, between —0.1a and —1.3c 
off the mean value. 

From (|58p one can see that smaller z/* (or ) allows smaller A resulting in higher 
values for \Jnl\- The thermal inflation scenario [26] after R 2 inflation may yield such 
values. 

Finally, let us comment on the SM Higgs inflation with nonminimal coupling to 
gravity [HI], which is also preferred by the Planck data [32] . As was emphasized in 
[35], predictions of n s and r from the Higgs inflation and R 2 inflation are almost 
degenerate. Only the slight difference comes from the reheating stage. In fact, the 
SUSY flat direction alters reheating significantly from that after the original R 2 
inflation, thereby producing non- negligible Jnl as opposed to the Higgs inflation. 
In our scenario the reheating temperature can be as high as 10 14 GeV and spatially 
modulated. 
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A Scalaron interaction with matter fields 



In generic /(i?)-type gravitational theory, 



S = A / f(R)^9d A x + S m , (64) 



the extra scalar degree of freedom, <p, of gravity is extracted by 

ip = ^M Pl ln\f'(R)\, (65) 



and its mass M s is given by M s = 1/ \/3f"(R). So in the original i? 2 -inflation model 
with f(R) = R + R 2 /(6M 2 ), the scalar on mass is constant and given by M which is 
fixed by the observed amplitude of curvature fluctuations. The supergravity model 
discussed in this paper has a more complicated action and the scalaron, which is in 
fact a complex field and we are focusing on its real part throughout [21] , has different 
masses at different regimes. 

Separating this scalaron degree of freedom through appropriate conformal rescal- 
ing, the rest of the gravitational action is given by the familiar Einstein term only. 

_KA£_ 

The action for a massive scalar field a, after the conformal rescaling a — > e v^a 
reads 



1 KCT K 2 CJ 2 

Scalar = --d^ad^a - -j=d ' m cj 'ip —d^ipd^ip 



^e-^V. (66) 

£4 



Similarly, the action for a canonical fermion is given by 



£fcrmion = -Ipplp ~ & ^^m^flp, (67) 

after rescaling ip — > e 4^^. The spinor covariant derivative is defined, for a U(l) 
gauge field A^, as 

prP = e» a7 a (d^ + T M - igA^, (68) 
r M = lh a ,^]e\y^, (69) 



Kip 



where is a tetrad (vierbein) field, rescaled as — > e ^e^ a , and is a 
spin connection for the Dirac spinor [TT]. For a U(l) gauge field A^ and a charged 
massless complex scalar <i> = xe^/v^, the transformed action is given by 



^vector = ~\f^F, u - \gTdnXdvX ~ *^=<T %X%V> ~ ^<T<V^ 



-yxV^K, (70) 
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where F^ u = d^V v - d v V^ = d^A u - d^A^, a gauge invariant field V /i = A li + <9 M C has 

Kip 

been denned, and rescaled as x — ► e ^X- If the gauge symmetry has been broken 
explicitly or spontaneously before rescaling, one could consider the gauge-breaking 
mass term 

w=4 e ^ y2 - (7i) 

The scalaron ip can decay into the scalar a, fermion and massive vector V~ M 
via trilinear interactions: 

where we have integrated by parts, used S(y/— <7>C sca i ar ) / 5a = 0, and left only linear 
terms in 93 from Eqs. ([66]) . (|67|) and ([7T|) . 

The decay rates of the scalaron are then given by [16j [T71 [11] 

r (^™) = mirMp \M. V'^TJ ' (73) 

2 \ 3/2 



5 



_> ^) = y » ' l £ , (74) 



where AC, AC, and Ay are the number of modes for each field. Note that we have 
added gauge-breaking mass term (jTTI) to the vector action (jTOj) ; otherwise the decay 
channel to vector modes do not show up classically. 

We are interested in the case masses of the decay products are given by a spatially- 
dependent expectation value of a scalar field <j), which is either the Higgs field in the 
SM or a flat direction field in SUSY, as m a = hep, = y4>, and my = gcj) with h, y, 
and g representing, scalar, Yukawa, and gauge couplings, respectively. Hereafter we 
consider the case masses of the decay products are much smaller than the parent 
particle, so that the phase space suppression factors in the above decay widths are 
negligible. Then we can classify each mode according to the dependence on <fi as 
follows. 

w , KM* , K(M>) 2 M a , AW) 4 

F[ip — > 00) = 2 -I 2 h 



192irM pl 48irM pl 4SnM pl M s 

= K(?SK + T S km + T 3M ), (76) 

- A^(#) 2 M S _ 
T(ip -> i/n(>) = ^ 87rM 2 — = M^Tfm, (77) 
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Here T$k from scalar kinetic term is independent of <j) and provides the main part 
of the decay width, while T$km and Tfm give equal contribution if h = y which is 
the case in SUSY models. The former is the interference contribution between scalar 
kinetic and mass terms, and the latter is due to the fermionic mass term. Finally 
Tsm and Tvm, both from their respective mass terms, are higher order in the mass 
ratio compared with Tskm and Tfm, so that their contribution will be minor, and 
we do not incorporate them to the subsequent analysis, which turns out to be a good 
approximation in our case with a large scalaron mass. Note that the equivalence of 
the above scalaron-decay approach with the Bogoliubov transformation approach to 
calculate particle production rate have been shown in [12\ [T4"] in the present context 
and in |31j in a different context. 
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